Equivalent synchronization of chaos in driven and in autonomous systems 
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It is shown that the synchronization behavior of a system of chaotic maps subject to either an 
external forcing or a coupling function of their internal variables can be inferred from the behavior 
of a single element in the system, which can be seen as a single drive-response map. From the 
conditions for stable synchronization in this single driven-map model with minimal ingredients, we 
find minimal conditions for the emergence of complete and generalized chaos synchronization in both 
driven and autonomous associated systems. Our results show that the presence of a common drive 
or a coupling function for all times is not indispensable for reaching synchronization in a system of 
chaotic oscillators, nor is the simultaneous sharing of a field, either external or endogenous, by all 
the elements. In the case of an autonomous system, the coupling function does not need to depend 
on all the internal variables for achieving synchronization and its functional form is not crucial for 
generalized synchronization. What becomes essential for reaching synchronization in an extended 
system is the sharing of some minimal information by its elements, on the average, over long times, 
independently of the nature (external or internal) of its source. 

PACS numbers: 05.45.-a, 05.45.Xt, 05.45.Ra 



I. INTRODUCTION 

Chaos synchronization is a fundamental phenomenon 
in dynamical systems (Pecora & Carroll 1990, Pikovsky 
et al. 2002). Its study has provided insights into many 
natural processes and motivation for practical applica- 
tions such as secure communications and control of dy- 
namical systems (Boccaletti et al. 2002, Uchida et al. 
2005, Argyris ei a/. 2005). Chaos synchronization is com- 
monly found in unidirectionally coupled systems, where a 
distinction can be made between a drive or forcing sub- 
system and another driven or response subsystem that 
possesses chaotic dynamics (Pikovsky et al. 2002). Com- 
plete synchronization occurs when the state variables of 
the drive and the response subsystems converge to a sin- 
gle orbit in phase space. More generally, generalized 
synchronization of chaos arises when a functional rela- 
tion different from the identity is established between the 
drive and the response subsystems (Rulkov et al 1995, 
Abarbanel et al. 1996, Kapitaniak et al. 1996, Hunt et 
al. 1997, Parlitz & Kokarev 1999, Zhou & Roy 2007). 

The auxiliary system approach (Abarbanel et al. 1996) 
shows that when a response and a replica subsystems are 
driven by the same signal, then the orbits in the phase 
spaces of the response and replica subsystems become 
identical and they can evolve on identical attractors, if 
their initial conditions lie on the same basin of attraction 
of the driven- response system. By extension, an ensemble 
of identical chaotic oscillators can also be synchronized by 
a common drive. In this case, generalized synchroniza- 
tion implies that the elements in the ensemble become 
synchronized into a single trajectory that evolves differ- 
ently from that of the drive. The specific functional form 
of the drive is not determinantal; the basic mechanism 



that leads to synchronization is the sharing of the same 
input by the oscillators for all times. 

On the other hand, there has been interest in the in- 
vestigation of chaotic synchronization and other forms 
of collective behaviors emerging in autonomous systems 
consisting of networks of mutually interacting nonlinear 
units without the presence of external influences (New- 
man et al. 2006, Manrubia et al. 2004). Synchroniza- 
tion in globally coupled oscillators is relevant in many 
chemical and biological systems and have been experi- 
mentally investigated (Wang et al. 2000, S. Dc Monte et 
al. (2007), Taylor et al. 2009). In autonomous systems 
where the units are coupled through their mean field, the 
units and the mean field can synchronize to a common 
chaotic trajectory (Kaneko 1989). This process can be re- 
garded as the analogous of complete synchronization in 
driven systems of elements subject to a common drive. 
Recently, the concept of generalized synchronization of 
chaos occurring in driven systems has been applied to 
the context of autonomous systems (Alvarcz-Llamoza & 
Cosenza 2008). This means that, under some circum- 
stances, the chaotic state variables in an autonomous 
dynamical system can be synchronized to each other but 
not to a shared coupling function containing information 
from those variables. 

The occurrence of both forms of chaos synchroniza- 
tion in driven and in autonomous systems suggests that 
the nature of the common influence acting on the ele- 
ments in a system is irrelevant; it could consist of an 
external forcing, or an autonomous interaction function 
(Alvarez-Llamoza & Cosenza 2008). At the local level, 
each element in the system is subject to a field, either ex- 
ternal or endogenous, that eventually induces some form 
of synchronization between that field and the element. 



Thus, the local dynamics can be seen as a single element 
subject to a drive evolving as the common field in each 
case. 

In this paper wc present a simple and unified scheme 
for the study of chaos synchronization in classes of driven 
and in autonomous dynamical systems. Our method is 
based on the reported analogy between a single driven 
map and a system of globally coupled maps (Parravano 
& Cosenza 1998). This analogy provided an explana- 
tion of dynamical clustering in globally coupled systems 
(Parravano & Cosenza 1998, Manrubia et al. 2004) and 
of stability of steady states in systems with delayed inter- 
actions (Masoller & Marti 2005). We show that the syn- 
chronization behavior of a system of chaotic maps sub- 
ject to either an external forcing or a coupling function of 
their internal variables can be inferred from the behavior 
of a single driven map in the system. When this single 
drive-response system reaches cither complete or gener- 
alized synchronization, the local analogy implies that an 
ensemble of identical maps subject to a common (exter- 
nal or internal) field that behaves as this drive, should 
also synchronize in a similar fashion. We show that the 
conditions for stable synchronization in a single driven 
map also describe the regions where synchronized states 
occur in both driven and autonomous extended systems. 

By studying a single drivcn-map model with minimal 
ingredients, we find minimal conditions for the emergence 
of chaos synchronization in both driven and autonomous 
associated systems. Our results show that the presence 
of a common drive or a coupling function for all times 
is not indispensable for reaching synchronization in an 
extended system of chaotic oscillators, nor is the simul- 
taneous sharing of the drive or the coupling function by 
all the elements in the system. In the case of an au- 
tonomous system, the coupling function does not need 
to depend on all the internal variables for achieving syn- 
chronization and, in particular, its functional form is not 
crucial for generalized synchronization. What seems es- 
sential for reaching synchronization in a system is the 
sharing of some minimal information by the elements, on 
the average, over long times, independently of its source. 

In Sec. II, the single driven map model is presented 
and the stability conditions for complete and generalized 
chaotic synchronization are analyzed. The synchroniza- 
tion behavior of a system of homogeneously driven maps 
is explored in Sec. Ill and that of a system of maps sub- 
ject to an asynchronous drive is shown in Sec. IV. An 
autonomous system of homogeneously coupled maps is 
considered in Sec. V. The behavior of an autonomous sys- 
tem of heterogeneously coupled maps is shown in Sec. VI. 
Conclusions are presented in Sec. VII. 



II. SINGLE DRIVEN MAP 

Our work is motivated by the practical aspect of 
searching for minimal requirements for the emergence of 
chaos synchronization in dynamical systems. As a mini- 



mal model, we consider a single map driven with a prob- 
ability p. 
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where f{xt) and g{yt) describe the dynamics of the driven 
and the drive variables, xt and yt, respectively, at discrete 
time t and the coupling relation between them is chosen 
to be of the diffusive form 



w{xt,yt) = (1 -e)f{xt) + eg{yt), 



(2) 



where e is the coupling strength. Since we are focused 
on chaotic synchronization, for the driven dynamics we 
shaU choose f{xt) = h+\n \xt\^ where & is a real parameter 
and St S (—00,00). This logarithmic map exhibits robust 
chaos, with no periodic windows and no separated chaotic 
bands, on the interval h £ [—1,1] (Kawabe & Kondo 
1991). Robustness is an important property in practical 
applications requiring reliable operation under chaos in 
the sense that the chaotic behavior cannot be modified 
by small perturbations of the system parameters. 

The linear stability condition for synchronization is 
determined by the Lyapunov exponents of the two- 
dimensional system Eq. ([1]). These are defined as 
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Lx and Ly arc the magnitude of the eigenvalues of 

[Ilt^o •^[xtTyt)Y^"^ , and J{xt,yt) is the Jacobian matrix 
for the system Eq.([T]), calculated along an orbit. A given 
orbit {xt, yt} from i^Otot^^T— 1 can be separated in 
two subsets, according to the source of the Xt variable, 
either coupled or uncoupled, that we respectively denote 
as A = {{xt,yt} : Xt = w{xt-i,yt-i)} possessing pT ele- 
ments, and B = {{xt,yt} ■ Xt = f{xt-i)} having {l-p)T 
elements. We get 

i/r 



/T-l 



l/T 



/ n ^- n /'(^*) ^ \ 






T-l 



V 



t:XteA t-.XtEB 

n 9'iyt) 



t=0 



(3) 

where Wx = ^ — (1 — e)/'(a:), and K is a poly- 
nomial whose terms contain products of Wx, e, and 
g'{yt) to be evaluated along time. Then Lx — 



xteA 



".n.,eB/'(-^*)]^/^ and Ly = [Ui^„'g'{yt)]'^ 



Thus we get 
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where Xg is the Lyapunov exponent of the map g{yt)- 
We shall assume chaotic drives g{yt)\ thus Ay > 0. Syn- 
chronization takes place when the Lyapunov exponent 



corresponding to the driven map is negative (Rulkov et 
al. 1995); i.e., Aj. < 0. For a given set of parameter val- 
ues, there is a definite value of the probability p at whieh 
the exponent A^^ may change its sign, from positive to 
negative, signaling the onset of synchronization and the 
appearance of a contracting direction in the dynamics of 
the single driven map system Eq. ([1]). Thus, the drive 
just needs to act a fraction p of the time in order to 
achieve synchronization in the system Eq. ([1]). 

When g = f, the condition A^; < implies complete 
synchronization, where xt ~ yt- In this case we get 
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A^=pln|l-e| + A/, 
where A/ is the Lyapunov exponent of the map /. 



(6) 
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FIG. 1: Boundaries Aa, = on the plane (p, e) for syn- 
chronization in the single driven map Eq. ([l| with / = 
— 0.7 -|- In Ixtj. (a) Complete synchronization (CS) with g = 
/ = —0.7 + ln\yt\. (b) Generalized synchronization (GS) with 

g{yt)^0.5 + ln\yt\^f. 

Figure [Ua) shows the stability boundary, given by 
Ax = 0, for the completely synchronized states of the 
system Eq. ^ on the space of parameters (p, e) for a 
drive g ^ f. On the other hand, if g ^ /, the condi- 
tion Aa; < corresponds to generalized synchronization, 
characterized by Xt ^ yt- Figure [IJb) shows the sta- 
bility boundary, given by A^, = 0, for the generalized 
synchronized states of the system Eq. ([T]) on the space of 
parameters (p, e) for a drive <? ^ /. 



III. HOMOGENEOUSLY DRIVEN SYSTEMS 

The auxiliary system approach (Abarbanel et al. 1996) 
implies that a driven map can synchronize on identical 
orbits with another, identically driven map. By exten- 
sion, let us consider a system of N intermittently and 
homogeneously driven maps, defined as 
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{l-e)f{x\)+eg{yt), 



with prob. p, 
with prob. (1 — p) 

(7) 
1,2,..., N) gives the state of the ith map 
at discrete time t, / describes the local dynamics, and e 
is the strength of the coupling to the drive g{yt)- 



where x\ {i 



Since each map in the system Eq. (O is subject to 
the same external influence (or lack of it) at any time, 
the properties of this system can be characterized from 
the behavior of the individual local dynamics. Thus, the 
occurrence of stable synchronization in the single driven 
map Eq. ([T} should lead to synchronization in the ex- 
tended system of maps Eq. ([T]) , even when the drive acts 
intermittently in both cases. A completely synchronized 
state in the system Eq. ([7]) is given by xl — xt = yt, 
\/i, and it can occur when g = f ■ On the other hand, if 
.9 7^ /j generalized synchronization, characterized by the 
condition x\ = xt ^ yt, Vi, may also arise in this system 
for p < 1. 

Synchronization in an extended system can be charac- 
terized by the asymptotic time-average (cr) (after discard- 
ing a number of transients) of the instantaneous standard 
deviations at of the distribution of state variables x\ , de- 
fined as 
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where St is the instantaneous mean field of the system. 



1 ^ 



(9) 



4=1 



Stable synchronization corresponds to (a) ~ 0. Here we 
use the numerical criterion (a) < 10"''. 

Figure [5] shows at as a function of time for the inter- 
mittently driven system Eq. (O subject to chaotic drives 
g ~ f and g ^ f , corresponding to complete and gener- 
alized synchronization, respectively. 
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FIG. 2: (Jt vs. t for the systems of driven maps Eq. ((Tj) with 
f{x\) ~ —0.7 + InliJI, for different forms of the drive g{y). 
(a) g — — 0.7 + lnjyt| ~ f', complete synchronization, (b) 
diyt) = 0.5-f ln|2;t| 7^ /; generalized synchronization. In 
both cases, A'^ = 10*, e — 0.6, p — 0.6, initial conditions were 
randomly and uniformly distributed such that Xi G [—10,3] 
and the drive g{yt) is applied starting at i = 1000. 



For a given value of the coupling strength e, there is 
a threshold value of the probability p required to reach 
either type of synchronization. In each case, the con- 
dition Aa; < for synchronization in the single driven 
map Eq. ([T]) implies stable complete or generalized syn- 
chronization in the intermittently driven system Eq. ([7]) . 
For the same values of parameters as in the single driven 
map, the region for the complete synchronized state on 
the space of parameters (p, e) of the system Eq. ([7]) with 
a drive g = f is the same as in Figure HJa) . Similarly, 
Figure [TJb) describes the region for generalized synchro- 
nization corresponding to the system Eq. ([7]) subject to 
a drive g ^ f- 



IV. HETEROGENEOUSLY DRIVEN SYSTEMS 

The equivalence between a single driven map and a sys- 
tem of driven similar maps also suggests that, under some 
circumstances, the collective behavior of an extended sys- 
tem of interacting elements can be inferred by consider- 
ing the dynamics of a single element at the local level. 
As an application of this idea, we consider a system of 
heterogeneously driven maps defined as 



'-t+i 



(1 - €.)f{xl) + eg{yt), with probability p, 
f{xl), with probability (1 — p). 

(10) 
The parameter p is the probability of interaction of a map 
with the drive g at a time t. The driven elements are 
randomly chosen with a probability p, so that not all the 
maps in the system receive the same external influence at 
all times. Thus, the average fraction of driven elements in 
the system Eq. (|10p at any given time is p. In comparison, 
the forcing of the elements in the homogeneously driven 
system Eq. ([7]) is simultaneous and uniform; each map 
receives the same influence from the drive g at any t 
with probability p. 

When the heterogeneously driven system Eq. ([T0| gets 
synchronized, we have xj = Xt- However, the synchro- 
nized solution exists only if g ~ f. Therefore, only com- 
plete synchronization xl = Xt — yt can occur in this 
system. Figure [3] shows at as a function of time for the 
system Eq. ([T0|) subject to a drive g = f. 

At the local level, each map in the heterogeneously 
driven system Eq. ([TU)) is subject, over long times, to an 
external forcing g with probability p. Thus, as in the 
homogeneously driven system Eq. ([7]), the synchroniza- 
tion behavior of system Eq. (jlOp can also be studied from 
the behavior of the single map driven with probability p, 
Eq. ([T]). In particular, if the system of maps Eq. ^TU\\ 
driven with g ~ f reaches a complete synchronized state 
for some values of parameters, then for this same set of 
parameters the single driven map Eq. ([T]) subject to the 
same drive should eventually exhibit a synchronized state 
similar to that of system Eq. ([T0| . Thus, the condition 
Aa; < for complete synchronization in the single driven 
map Eq. ([T]), that implies stable complete synchroniza- 
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FIG. 3: at vs. t for the systems of driven maps Eq. (|10p 
with f(xl) = -0.7 + ln|a;j|, for g = -0.7 + ln|yt|; complete 
synchronization. Parameters are A'^ — W^, e = 0.6, p = 0.6; 
initial conditions were randomly and uniformly distributed 
such that Xi G [—10, 3] and the drive g{yt) is applied starting 
at t = 1000. 



tion in the homogeneously driven system Eq. ([7]), also 
describes the stability of the complete synchronized state 
xl = Xt — yt, yi, in the heterogeneously driven system 
Eq. (Uni). The stability boundary A^^ = in Fig. [Ha) 
for the driven map with g = f coincides with both, the 
boundary that separates the region where complete syn- 
chronization occurs on the space of parameters (p, e) for 
the homogeneously driven system Eq. ([7]) and the bound- 
ary for complete synchronization in the heterogeneously 
driven system Eq. ([TUl) . 



V. AUTONOMOUS SYSTEMS WITH 
HOMOGENEOUS COUPLING 

The dynamics of the single driven map Eq. ^ can 
also be compared with the dynamics of an autonomous 
system of maps that share a global interaction with a 
probability p, 
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(1 - e)fixl) + eHixl : j E Qt), prob. p, 



prob. (1 — p) 

(11) 
where e represents the magnitude of the coupling to the 
global interaction function H; and Qt is a subset hav- 
ing q < N elements of the system that may be chosen 
at random at each time t. Each map receives the same 
information from the feedback coupling function H at 
any t with probability p. When this autonomous sys- 
tem gets synchronized at some values of parameters, we 
have xl = Xt- Thus for this same set of parameters, 
the single driven map subject to a forcing that satisfies 
giut) = ii(x\ = Xt '■ j G Qt) for long times should exhibit 
a synchronized state similar to that of the autonomous 
system Eq. (jlip . Thus, besides complete synchronization 
where xl = Xt = H , other synchronized states, character- 
ized by xl = Xt j^ H should also occur in the autonomous 
system, Eq. (jlip . for appropriate values of parameters. 
These states can be described as generalized synchroniza- 
tion in autonomous dynamical systems (Alvarez-Llamoza 
& Cosenza 2008). 

As an example of complete synchronization in the sys- 



tern Eq. ((TT|) , consider a partial mean field coupling func- 
tion defined as 



H{xi:jEQt)^-J2fi^i)^ 



(12) 
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where q < N maps arc randomly chosen at each time 
t. In this case the autonomous system Eq. pT|) can be 
expressed in vector form as 



Xt+l 



(l-e)I+-Gt 



f (xt), with probability p, 
I f(xt), with probability (1 —p), 



(13) 
where the A^-dimcnsional vectors xt and f (xt) have com- 



ponents [Xt]i 



and [f(xt)].i = /(xj), respectively, I is 



the N X N identity matrix, and Gt is an TV x A^ matrix 
that at each time t possesses q randomly chosen columns 
that have all their components equal to 1 while the re- 
maining N — q columns have all their components equal 
to 0. The case q = N and p = 1 corresponds to the 
usual mean field global couphng (Kaneko 1989). The 
linear stability analysis (Waller & Kapral 1884) of the 
complete synchronized state f{xl) = ,f{xt) — H yields 



(l-e) + -afc 

q 



e^/ 



< 1. 



(14) 



where a^ = ^ofc? {k = 0, 1 . . . , A^ — 1) are the set of eigen- 
values of the matrix Gt for any f , with the zero eigenvalue 
having [N — l)-fold degeneracy. The eigenvector corre- 
sponding to fc = is homogeneous. Thus only perturba- 
tions of Xt along the other eigenvectors may destroy the 
coherence. Thus, condition Eq. p^ with k ^Q becomes 



pln|l-e| + A/ <0, 



(15) 



which is the same condition for stability of complete syn- 
chronized states in the driven map, Eq. (|6]), when g = f. 
Thus the boundary that separates the region where com- 
plete synchronization occurs on the space of parameters 
{p, e) for the autonomous system Eq. (|lip , with H given 
by Eq. ([T2|) for any value of g, coincides with the stabil- 
ity boundary A^; = in Fig. [5] for the driven system with 

9 = f- 

Figure mja) shows the bifurcation diagram of the mean 
field St as well as (a) as functions of the probability p 
for the autonomous system with the coupling function H 
given by Eq p^. with a fixed value of e. The mean field 
St is chaotic for all values of p. There is a critical value of 
p above which (a) = and /(xj) = f{xt) = H, indicating 
that the system is completely synchronized in a chaotic 
state for that range of the probability of interactions. 

For other functional forms of the coupling function H 
it is possible to find generalized synchronized states in 
the autonomous system Eq. (jlip . For example, consider 
the coupling function 
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FIG. 4: Left vertical axis: bifurcation diagram of St as a 
function of p for the homogeneously coupled autonomous sys- 
tem Eq. (|11|) . Right vertical axis: (p-) vs p, continuous line. 
For each value of p, St and {a) were calculated at each time 
step during a run starting from random initial conditions on 
the local maps, uniformly distributed on the interval [—10, 3], 
after discarding the transients, (a) Ht given by Eq. (|12|) 
with q = 0.5A^; the region where the complete synchroniza- 
tion state occurs is labeled CS. (b) Ht given by Eq. (|16|) 
with k = 1.2 and q — 0.5N; the region where generalized 
synchronization occurs is indicated by GS. In both cases, 
f(x) = -0.7 + log \x\, e^ 0.5, and N = 10**. 



where A: is a constant and q < N elements are chosen at 
randonr at each time t. 

Figure HJ^b) shows both the bifurcation diagram of St 
and the quantity (a) as functions of the probability p for 
the autonomous system Eq. (fTTj) with the coupling func- 
tion H given by Eq P^ . with a fixed value of e. Again, 
there is a value of p above which (a) vanishes, signaling 
the onset of a synchronized chaotic state. In this region 
f{xl) — f{xt) ^ H = k + f{xt), independently of the 
value of q, corresponding to generalized synchronization 
in the autonomous systenr Eq. (|11[) . 

We have verified that the regions of the parameter 
space {p, e) where complete or generalized synchronized 
states emerge in the autonomous system with intermit- 
tent global interactions, Eq. (fTTj). are the same regions 
that correspond to those states for the single driven nrap 
in Figure [TJ 



VI. AUTONOMOUS SYSTEMS WITH 
HETEROGENEOUS COUPLING 

In analogy to a system of heterogcncously driven maps, 
Eq. (jlOp . we can consider the following autonomous cou- 
pled map system, 



^t+i 



{l-e)f{xl) + eH{xi:jeQt), prob. p, 
f{xl), prob. (1-p). 

(17) 
The parameter p is the probability of interaction of an 
element with the coupling function at time t. Thus the 
average fraction of fully connected elements at any given 
time is p. For p = 1, the autonomous systems Eq. (J17p 
and Eq. ([TT|) are identical. 

Not all the maps in the system Eq. pT)) experience 
the same global interaction at a given time. However, at 
the local level, each map in this autonomous system is 
subject, on the average, to the same coupling function 
H with probability p over long times. Similarly, each 
map in the autonomous homogeneously coupled system 
Eq. (|lip experiences the coupling H with a probability 
p in time. This is also analogous to the condition of the 
local dynamics in the heterogeneously driven system of 
maps, Eq. (jlO[) . if g exhibits the same temporal evolution 
as H. Then, the synchronization behavior of the system 
Eq. ((T7)) should also be similar to the behavior of the 
single map driven, Eq. ([T]) subject to a drive that satisfies 
g{yt) = H for asymptotic times. 

When the autonomous system Eq. pT]) reaches a syn- 
chronized state, we have f{x\) = f[xt)- However, this 
synchronized solution exists only if _ff = f{xt). There- 
fore, as in the case of the heterogeneously driven system, 
only complete synchronization {x\) = f{xt) = H arises 
in the system Eq. ((T7)) . 

As an example of complete synchronization in the au- 
tonomous system Eq. ([T7)) . we choose the form Eq. p^ 
as the coupling function H{xl : j £ Qt) Figure [5] shows 
St and the quantity (ct) as functions of the probability p, 
for a fixed value of e. For the given values of parameters, 
there is a threshold value of p at which {a) = 0, indi- 
cating that the elements become synchronized; in this 
case, complete synchronization, f{x\) = f{xt) ~ H, in- 
dependently of the value of q. Thus, the reinjection of 
a coupling function containing partial information of the 
system to a fraction of randomly selected elements suf- 
fices to achieve complete synchronization in the system. 

Beyond the common regions for complete or gener- 
alized synchronization, the driven and the autonomous 
systems may exhibit different collective behaviors in the 
parameter space. This is particularly manifested in the 
autonomous system with inhomogeneous coupling. Be- 
sides complete synchronization, the instantaneous mean 
field St in Figure [5] (left vertical axis) shows more complex 
collective behaviors arising in the system Eq. (|17|) . For 
the parameter value b = —0.7, the local maps are chaotic. 
However, for some range of the parameters p and e, St 
reveals the existence of global periodic attractors. In Fig- 
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FIG. 5: Left vertical axis: bifurcation diagram of St as a 
function of p for the partially coupled autonomous system 
Eq. dlTl) with Ht defined in Eq. dH} and q = 0.57V. The region 
where the complete synchronization state occur is labeled CS. 
Right vertical axis: (a) vs p, continuous line. For each value 
of p, St and (a) were calculated at each time step during a 
run starting from random initial conditions on the local maps, 
uniformly distributed on the interval [—10, 3], after discarding 
the transients. Here, f{x) = — 0.7 + log |a;|, e = 0.4, and 
iV = 10^ 



urc El collective periodic states of St at a given value of 
the parameter p appear as sets of small vertical segments 
which correspond to intrinsic fluctuations of the periodic 
orbits of the mean field. Global attractors of period 1, 
2, 4, and 8 are observed in the mean field as the param- 
eter p is varied, and they become better defined as the 
system size N is increased. The amplitudes of the global 
periodic motions of St do not decrease with an increase 
of A^. As a consequence, the variance of St itself does 
not decay as A^^ with increasing A^ but rather it satu- 
rates at some constant value related to the amplitude of 
the collective period. This is a phenomenon of nontriv- 
ial collective behavior, where macroscopic quantities in 
a spatiotemporal dynamical system exhibit spontaneous 
ordered evolution coexisting with local chaos (Chate & 
Manneville 1998). This phenomenon arises in the het- 
erogeneously coupled autonomous system, rather than in 
the homogeneously coupled one. In the first case, the 
coupling function presents spatial fluctuations than can 
be interpreted as a distributed noise. In fact, it has been 
shown that local noise can enhance the emergence of or- 
dered collective behavior in coupled chaotic map systems 
(Shibata et al. 1999). Note that the emergence of col- 
lective periodic behavior in our case shown in Figure [5] 
cannot be attributed to the presence of periodic windows 
in the local dynamics since the logarithmic map possesses 
robust chaos for 6 G [1, 1] (Alvarcz-Llamoza et al. 2007). 
Although complete or generalized synchronization of 
chaos in both, driven and autonomous systems, can be 
characterized from the knowledge of the dynamical re- 
sponse of a single driven map, the transient behavior to 
reach such a state may be different in each case. Fig- 
ure ^ shows the average time Tg required to attain com- 
plete or generalized synchronization as a function of p in 



the driven and the autonomous systems considered in 
this paper. The time Tg for generahzed synchronization 
is larger than that for complete synchronization in an 
extended system subject to a field, either external or en- 
dogenous. In general, a homogeneous, synchronous field 
acting on a system is more efficient for achieving complete 
synchronization. The curves Tg vs. p for complete syn- 
chronization in either a driven or an autonomous system 
subject to a homogeneous field are practically indistin- 
guishable. The same curve is also obtained for a single 
driven map possesing similar parameter values, as one 
may expect from the local dynamics analogy. Varying 
the size of the system N does not significantly affect the 
curves in Figure ([6|). 
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FIG. 6: Average time Ts to reach synchronization as a func- 
tion of the probability p. The error bars on each curve corre- 
spond to standard deviations resulting from 100 realizations 
of random initial conditions for each value of p. (a) Driven 
systems with fixed e = 0.4 and A'^ = 10 . Top curve: ho- 
mogeneously driven system, Eq. (O, with / = —0.7 -I- log |a;|, 
g = 0.5 + log \y\ (generalized synchronization). Middle curve: 
heterogeneously driven system, Eq. (|10|l . with f = g = —0.7-1- 
log|a;| (complete synchronization). Bottom curve: homoge- 
neously driven system, Eq. ((Tj, with f = g = —0.7 -I- log \x\ 
(complete synchronization), (b) Autonomous systems with 
fixed e = 0.4 and N = 10*. Top curve: autonomous homoge- 
neously coupled system, Eq. Hi}, with H given by Eq (|16|l . 
q = 0.5A and k — 1.2 (generahzed synchronization). Middle 
curve: autonomous heterogeneously coupled system, Eq. (|17|l , 
with H given by Eq. H12|) . q — 0.5N (complete synchroniza- 
tion). Bottom curve: autonomous homogeneously coupled 
system, Eq. (fTTj) . with H given by Eq. (fT2)) . q = 0.5N (com- 
plete synchronization). 

Note that for achieving synchronization in cither the 
homogeneous or heterogeneous autonomous system, the 
coupling function H does not need to depend on all the 
internal variables; nor must it be active for all times or 
act simultaneously on all the elements. What becomes 
essential for synchronization is that all the elements in an 
autonomous system share some minimal information, on 
the average, over time. This is the same process that 
leads to synchronization in the associated driven sys- 
tems Eq. ([7]) and Eq. ([TU)) . respectively, and it is also 
the mechanism for synchronization in the single driven 
map, Eq. ([T}. Thus, the nature of the common infiuence, 
either external or internal, affecting the local units be- 
comes irrelevant for the emergence of synchronization in 



driven or autonomous systems. 



VII. CONCLUSIONS 

We have shown that the synchronization behavior of 
system of chaotic maps subject to cither an external forc- 
ing or a coupling function of their internal variables can 
be inferred from the behavior of a single element in the 
system. The local dynamics can be seen as a single driven 
map. When this single drive-response system reaches ei- 
ther complete or generalized synchronization, the local 
analogy implies that an ensemble of identical maps sub- 
ject to a common field, cither external or internal, that 
behaves as this drive, should also synchronize in a similar 
manner. 

By considering the conditions for stable synchroniza- 
tion in a single driven-map model with minimal ingredi- 
ents, we have found some minimal, equivalent conditions 
for the emergence of complete and generalized chaos syn- 
chronization in both driven and autonomous associated 
systems. 

For the homogeneously driven or coupled autonomous 
systems, the sharing of the infiuence by the elements 
takes place only a fraction of the time. In the case of 
the heterogeneously driven or coupled autonomous sys- 
tem, the influence is shared only by a fraction of ele- 
ments at any time. Our results show that the presence 
of a common drive or a coupling function for all times 
is not indispensable for reaching synchronization in an 
extended system of chaotic oscillators, nor is the simul- 
taneous sharing of the drive or the coupling function by 
all the elements in the system. In the case of an au- 
tonomous system, the coupling function does not need 
to depend on all the internal variables for achieving syn- 
chronization and, in particular, its functional form is not 
crucial for generalized synchronization. What matters 
for reaching synchronization in a system is the sharing 
of some minimal information by the elements, on the av- 
erage, over long times, independently of its source. This 
result suggests that these systems possess an ergodicity 
property. 

Future extensions of this work include the investigation 
of some quantity, such as the transfer entropy (Schreiber, 
2000) for measuring this minimal amount of information 
required for either type of synchronization, and the possi- 
bility of prediction of other forms of collective behaviors 
or patterns observed in dynamical networks, from the 
knowledge of the occurrence of these phenomena in an 
associated single driven-response system. 
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